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Summary 

A formal procedure for the probabilistic design evaluation of a composite structure is described. The 
uncertainties in all aspects of a composite structure (constituent material properties, fabrication variables, 
structural geometry, and service environments, etc.), which result in the uncertain behavior in the 
composite structural responses, are included in the evaluation. The probabilistic evaluation consists of: 
(1) design criteria, (2) modeling of composite structures and uncertainties, (3) simulation methods, and 
(4) the decision-making process. A sample case is presented to illustrate the formal procedure and to 
demonstrate that composite structural designs can be probabilistically evaluated with accuracy and 
efficiency. 


Introduction 

Composite materials are widely used in modem structures for high performance and reliability. 
However, because these structures usually operate in hostile and random service environments, it is 
difficult to predict the structural performance. In addition, experiments show that the composite structural 
behavior exhibits wide scatter as a result of the inherent uncertainties in design variables. The design 
variables, known as primitive variables, include the fiber and matrix material properties at the constituent 
level; fiber and void volume ratios, ply misalignment and ply thickness’ the fabrication process, and 
random structure size, boundary conditions, loadings, and environment. 

The scatter in the structural behavior cannot be computationally simulated by the traditional 
deterministic methods that use a safety factor to account for uncertain structural behavior; thus, the 
structural reliability cannot be discerned. A probabilistic design methodology is needed to accurately 
determine the structural reliability of a composite structure. 

NASA Glenn Research Center has developed a formal methodology to efficiently and accurately 
quantify the scatter in the composite structural response and to evaluate the composite structural design, 
while accounting for the uncertainties at all composite levels (constituent, ply, laminate, and structure) 
(fig. 1). This methodology, which integrates micro- and macrocomposite mechanics and laminate 
theories, finite element methods, and probability algorithms, was implemented through the computer code 
IPACS (Integrated Probabilistic Assessment of Composite Structures) (fig. 2) (ref. 1). IPACS is used to 
evaluate composite structures probabilistically for all types of structural performances such as instability, 
clearance, damage initiation, delamination, microbuckling, fiber crushing, and resonance damage. Since 
IPACS uses a special probability algorithm FPI (Fast Probability Integrator) (ref. 2) instead of the conventional 
Monte Carlo simulation, tremendous computational time can be saved. Therefore, a probabilistic composite 
structural analysis, which cannot be done traditionally, becomes desirable especially for large structures with 
many uncertain variables. A typical case is analyzed herein to demonstrate the code IPACS for the 
probabilistic evaluation of composite structures and to illustrate the formal design evaluation methodology. 
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Figure 1. — Concept of probabilistic assessment of composite structures. 



Figure 2. — Probabilistic design assessment of composite structures. 


Fundamental Approach 

The fundamental approach for the probabilistic evaluation is as follows: 

(1) Identify all possible uncertain variables at all composite scale levels. 

(2) Determine the probabilistic distribution function (PDF) for each variable. 

(3) Process all random variables through an analyzer that consists of micro and macrocomposite 
mechanics, laminate theories, structural mechanics, and probability theories. 

(4) Extract useful information from the output of the analyzer and check against the probabilistic design 
criteria. 

The uncertainties in a composite structural design can originate at different composite scale levels. At 
the constituent level, the material properties for the fiber and matrix are the major sources of 
uncertainties. Typical values are listed in table I. At all stages of the fabrication process the fabrication 
variables, such as fiber volume ratio, void volume ratio, ply misalignment and ply thickness, show 
considerable scatter. At the structure level, variation of the geometry during the assembly stage, uncertain 
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boundary conditions, and random thermal-mechanical loads contribute significantly to the scatter in the 
composite structural response. 


TABLE I.— MATERIAL PROPERTIES AT THE CONSTITUENT LEVEL 
FOR SKIN AND STRINGERS 


Property 


E fll , Mpsi 
E 02 , Mpsi 
Gfi 2 , Mpsi 
G 03 , Mpsi 

Vfl2 

V C3 

a f n, ppm/ °F 
a£ 22 , ppm / °F 

p f , lb/in . 3 
N f 

d f , in. 

C f , Btu/in. °F 
K fl i, Btu- in./hr in . 2 
°F 

K Q2 , Btu- in./hr in . 2 
°F 

Kgj, Btu- in./hr in . 2 
°F 

Sfx, ksi 
Sfc. ksi 
E m , Mpsi 
G m , Mpsi 
v m 

a m , ppm/ °F 
p m , lb/in . 3 
C m , Btu/in. °F 
K m , Btu- in./hr in . 2 
°F 

S mT , ksi 
S m c. ksi 
S mS , ksi 
p m ,(in./in.)/l% 
moist 

D m , in. 3 /sec 


Assumed distribution type 


Mean 3 


Assumed uncertainty scatter 


Nonnal 


Constant 

Nonnal 


Weibull 

Weibull 

Nonnal 


Weibull 

Weibull 

Weibull 

Nonnal 

Nonnal 


31.0 
2.0 
2.0 
1.0 
0.2 
0.25 
-0.55 
5.6 
0.063 
10 000 
0.0003 
0.17 
580 
58 
58 
400 
400 
0.5 
0.185 
0.35 
42.8 
0.0443 
0.25 
1.25 
15 
35 
13 
0.004 
0.002 


±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 
±5% from 


mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 


■“Typical values for graphite- fiber/epoxy-matrix composites at 0.6 fiber volume ratio. 


Probabilistic Composite Structural Design Evaluation 

This paper describes a probabilistic design evaluation methodology for composite structures with an 
acceptable or preassigned reliability. The details of the evaluation are described in the following three 
steps. 

Step 1: Probabilistic Design Criteria Setup: A typical design criterion can be stated as follows: The 

probability of a failure event should be less than an acceptable value, say 10' 3 . A failure event occurs 
when a structural response is greater than the allowable response. This probability is defined as the failure 
probability. The allowable response divides the possible response domain into safe and failure regions as 
shown in figure 3. The predicted failure probability is the area under the probability density function in 
the failure region. The critical response (fig. 3) is determined by IPACS such that the probability of a 
response exceeding this critical value is in the safe region. When the critical response falls within the safe 
region, the design is acceptable. When the critical response falls within the failure region, the design is 
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Probabilistic design criterion 
Minimum allowable performance X 
Probability (operating performance < X) < 0.001 

Critical performance Y 

Probability (operating performance < Y) = 0.001 


Allowable 



Performance 

Figure 3. — Safe and failure regions in a probability space. 

unacceptable and requires a redesign. Sample probabilistic design criteria for the various failure modes 
are described as follows: (1) Instability — The probability that the buckling load is smaller than the design 
load should be less than 10" 3 . (2) Clearance — The probability that the nodal displacement is greater than 
the allowable tolerance should be less than 10“ 3 . (3) Resonance avoidance — The probability that the 
natural frequency is greater that its upper bound should be less than 1 0’ 3 . (4) Delamination — The 
probability of delamination occurrence should be less than 10“ 3 . 

Step 2: Probabilistic Simulation Using IPACS: IPACS integrates several NASA in-house computer 

codes developed in recent years such as COBSTRAN (Composite Blade Structural Analyzer), PICAN 
(Probabilistic Integrated Composite Analyzer) and MHOST (Marc Hot Section Technology). 
COBSTRAN is a dedicated finite element model generator for structural analysis of composite structures. 
PICAN uses ICAN (Integrated Composite Analyzer) computer code for composite mechanics. This code 
has evolved over the last 40 years and has been verified with experimental data for all aspects of 
composites. PICAN enables the computation of the perturbed and probabilistic composite material 
properties at the ply and laminate levels. MHOST performs structural analyses using verified finite 
element methods. These analyses determine the perturbed and probabilistic structural response at global, 
laminate, and ply levels. PICAN and MHOST share the FPI module for the application of the fast 
probability integration algorithm to obtain cumulative distribution functions (CDF’s) of the material 
properties and the structural responses. 

FPI is an approximate technique for the probabilistic analysis of the structural performance and has 
the major advantage of speed. FPI techniques are orders of magnitude more efficient than Monte Carlo 
simulation methods. This is especially true in the tails regions of the distribution; that is, at very high or 
low probabilities since the FPI solution time is independent of the probability level. Conversely, in Monte 
Carlo simulation methods, the computational time increases with very high or low probability levels. 

Also, FPI allows evaluation of information that describes the relative importance of each random variable. 
These sensitivity factors can be a valuable aid in optimizing a design. 

In IPACS, the probabilistic evaluation of composite structures starts with the identification of 
uncertain primitive variables at constituent and ply levels. These variables are then selectively perturbed 
several times to create a data base. The data base is used to establish the relationship between the desired 
structural response (or the desired material property) and the primitive variables. For every given 
perturbed variable, micromechanics is applied to determine the corresponding perturbed mechanical 
properties at the ply and laminate levels. Laminate theory is then used to ascertain the resultant force- 
moment strain-curvature relationship. With this relationship at the laminate level, a finite element 
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perturbation analysis is performed to find the structural response that corresponds to the selectively 
perturbed primitive variables. This process is repeated until enough data are generated and the proper 
relationship between structural response and primitive variables can be established through a numerical 
procedure. 

Step 3: Decision Making and Redesign: IPACS simulates the PDF of a given structural response, such 

as buckling load, displacement, local stress, local strength, vibration frequencies, and fatigue life. The 
probability of a design violation for each criterion can be calculated with these PDF functions. When the 
failure probability is greater than the acceptable value, say 1CT 3 , the composite structural design should be 
rejected. To redesign a composite structure, one can use the sensitivity factors from the IPACS analysis. 
Sensitivity factors ra nk the random variables based on their contribution to this failure probability. 
Therefore, a redesign will be guided by this information with manufacturing control of the mean and the 
standard deviation (stdv) of the appropriate random variables. 


Demonstration Case and Discussion 

A stiffened composite cylindrical pipe is probabilistically designed against probabilistic design 
criteria. The cylindrical pipe is 2 ft. in diameter and 20 ft. long. The structure is modeled by 588 four- 
noded shell elements and 600 active nodes (6 degrees of freedom per node). The composite pipe consists 
of the skin, three horizontal circumferential frames, and four vertical stringers. The laminate 
configurations for the skin, frames, and stringers are [+45/02/±45/02/±45/0/90] s , [O 24 ] and [O 24 ], 
respectively. The pipe is assumed to be supported at one end by a set of translational and torsional spring 
constants and free at the other end. When the spring constant approaches infinity, a completely fixed 
boundary condition is simulated. When the spring constants are set to zero, a free boundary condition is 
simulated. For a given set of spring constants, a partially fixed boundary condition is modeled. The pipe 
is subjected to axial (F x ) and lateral (F y ) loads as well as torsional moments (M xx ) at its free end. 

The uncertain variables are identified at the constituent, ply, and structure levels. At the constituent 
level, 17 material properties for the graphite fiber and 12 for the epoxy matrix of the skin, frames, and 
stringers are assumed to be uncertain variables. Their probability distribution types and associated 
parameters are listed in table I. At the ply level, the fabrication variables (fiber and void volume ratio, ply 
orientation, and ply thickness) are treated as random variables. Their statistics are shown in table II. At 
the structure level, spring constants that simulate a partially fixed boundary condition are assigned by a 
probability distribution as are the loading conditions. Their statistics are shown in table III. The composite 
pipe was designed against two design criteria: clearance and delamination. The results showed that the 
uncertainty range in the end displacement is between 1 and 3 percent, of the pipe length and it is most 
sensitive to uncertainties in the wall thickness perimeter variables. The probable delamination failure is 
most sensitive to sheer strength of the wall thinness. 


TABLE II.— FABRICATION VARIABLES AT PLY LEVEL 


Variable 

Assumed distribution type 

Mean 

Assumed uncertainty scatter 

Fvr 

Normal 

0.60 

0.05 

wr 


0.02 

0.05 

0 P , deg 


0.00 

0.9 (stvd) 

tpsk, in. 


0.005 

0.05 

W> in. 


0.02 

0.05 


TABLE III.— UNCERTAINTIES IN STRUCTURAL LEVEL 


Uncertainty 

Assumed distribution type 

Mean 

Assumed uncertainty scatter 

K cT r, lb/in. 

Nonnal 

30xl0 6 

0.20 

K cT o, lb-in./rad 


12xl0 2 

0.20 

F x , kip 


288 

0.05 

F y , kip 


5.76 

0.05 

M xx , kip- ft 


576 

0.05 
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Concluding Remarks 

A formal methodology is described in this report for the probabilistic design evaluation of composite 
structures. This methodology, integrating micro and macrocomposite theory, structural mechanics (finite 
element methods), and probability algorithms, performs a probabilistic evaluation of composite structural 
designs considering all identifiable uncertain variables at all composite levels. Composite structural 
designs can be evaluated against specific probabilistic design criteria demonstrating that such designs can 
be computationally evaluated by using the probabilistic computer code IPACS. Information for an 
efficient design can also be obtained. Specifically, for the demonstration case, the uncertainty range in the 
end displacement was between 1 and 3 percent of the pipe length and was most sensitive to the 
uncertainties in the skin-related primitive variables. Conversely, probable delamination failure was most 
sensitive to the shear strength of the skin. 
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